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In pp scattering at LHC energies, large numbers of elementary scatterings will contribute signif- 
icantly, and the corresponding high multiplicity events will be of particular interest. Elementary 
scatterings are parton ladders, identified with color flux-tubes. In high multiplicity events, many of 
these flux tubes are produced in the same space region, creating high energy densities. We argue 
that there are good reasons to employ the successful procedure used for heavy ion collisions: matter 
is assumed to thermalizes quickly, such that the energy from the flux-tubes can be taken as initial 
condition for a hydrodynamic expansion. This scenario gets spectacular support from very recent 
results on Bose-Einstein correlations in pp scattering at 900 GeV at LHC. 



I. INTRODUCTION 

After one decade of RHIC experiments it seems to be 
certain that heavy ion collisions at RHIC energies pro- 
duce a new state of matter which expands as an almost 
ideal fluid [1-8], whereas proton-proton scattering is usu- 
ally considered to be a reference system, theoretically well 
under control via perturbative techniques. Although at 
very high energy, hadrons experience multiple scatterings 
when they hit protons or neutrons, inclusive cross sections 
calculations becomes quite simple due to the fact that 
different multiple scattering contributions cancel due to 
destructive interference (AGK cancellations). The cor- 
responding formulas are simple and can be expressed in 
terms of parton distributions functions, based on evolu- 
tions equations. 

However, in particular at LHC energies where we ex- 
pect large numbers of scatterings to contribute signifi- 
cantly, it becomes interesting to study event classes cor- 
responding to a large number of scatterings (in practice: 
high multiplicity events). Here, one needs partial cross 
sections, corresponding to a particular multiple scattering 
type (single, or double, or triple...). Gribov-Regge theory 
provides a solution, in particular when energy sharing is 
properly taken into account, as in the EPOS approach. 

High multiplicity events are very interesting for the fol- 
lowing reasons: in EPOS for example, a single scattering 
amounts to the exchange of a complete parton ladder, in- 
cluding initial state radiation. The whole object is identi- 
fied as a pair of color flux tubes, which finally break into 
many pieces (hadrons). In high multiplicity events, with 
many scatterings involved, we have many partons lad- 
ders participating, and therefore a large number of flux 
tubes sitting essentially on top of each other - as in heavy 
ion scattering at RHIC. In the heavy ion case, we simply 
compute the energy density corresponding to these flux 
tubes (from string theory), assume thermalization, and 



then perform a hydrodynamic expansion based on these 
initial conditions 0. 

Since the energy densities reached in high multiplic- 
ity proton-proton collisions are comparable to the ones 
achieved in gold-gold scattering at RHIC, we will apply 
the same procedure. The usual argument against this 
approach is the small size of the pp system, but since we 
know by now that the size of the space fluctuations in 
an even-by-event treatment in AuAu scattering is of the 
order of 1-2 fm, and AuAu seems to be driven by hydro- 
dynamic flow, there is no reason not to do so for high 
multiplicity pp. 

In this paper, we will briefly review the flux-tube/hydro 
approach of with special emphasis on pp scattering. 
After some elementary checks concerning particle distri- 
butions, we come to the main result of this paper: the 
hydrodynamic expansion modifies drastically the space- 
time behavior of the evolution, compared to basic pic- 
ture where the flux-tubes decay independently. And this 
space-time structure can be clearly "seen" when investi- 
gating Bose-Einstein correlations, and the recently pub- 
lished results from ALICE confirm the "hydrodynamic 
scenario". 



II. MULTIPLE SCATTERING 

A. Parton evolution 

An elementary scattering within the EPOS approach 
Q is given by a so-called "parton ladder", see fig. [TJ 
representing parton evolutions from the projectile and 
the target side towards the center (small x). The evo- 
lution is governed by an evolution equation, in the sim- 
plest case according to DGLAP. In the following we will 
refer to these partons as "ladder partons", to be distin- 
guished from "spectator partons". Such a parton ladder 
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Figure 1: (Color online) Elementary interaction in the EPOS 
model. 



may be considered as a longitudinal color field or flux- 
tube, conveniently treated as a relativistic string. The 
intermediate gluons are treated as kink singularities in 
the language of relativistic strings. This flux tube ap- 
proach is just a continuation of 30 years of very successful 
applications of the string picture to particle production 
in collisions of high energy particles (9t-[l2|. in particu- 
lar in connection with the parton model. An important 
issue at high energies is the appearance of so-called non- 
linear effects, which means that the simple linear par- 
ton evolution is no longer valid, gluon ladders may fuse 
or split. More recently, a classical treatment has been 
proposed, called Color Glass Condensate (CGC), having 
the advantage that the framework can be derived from 
first principles (l3l4T7j . Comparing a conventional string 
model like EPOS and the CGC picture: they describe 
the same physics, although the technical implementation 
is of course different. All realistic string model implemen- 
tations have nowadays to deal with screening and satura- 
tion, and EPOS is not an exception, see Without 
screening, proton-proton cross sections and multiplicities 
will explode at high energies. 

A phenomenological treatment of non-linear effects in 
EPOS employs two contributions: a simple elastic rescat- 
tering of a ladder parton on a projectile or target nucleon 
(elastic ladder splitting), or an inelastic rescattering (in- 
elastic ladder splitting), see fig. [5J The elastic process 
provides screening, therefore a reduction of total and in- 
elastic cross sections. The importance of this effect should 
first increase with mass number (in case of nuclei be- 



(a) nucleons 
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Figure 2: (Color online) (a) Elastic "rescattering" of a ladder 
parton. We refer to elastic parton ladder splitting; (b) In- 
elastic "rescattering" of a ladder parton. We refer to inelastic 
parton ladder splitting. 

ing involved), but finally saturate. The inelastic process 
will affect particle production. Both, elastic and inelas- 
tic rescattering must be taken into account in order to 
obtain a realistic picture. 

To include the effects of elastic rescattering, we first 
parametrize a parton ladder (to be more precise: the 
imaginary part of the corresponding amplitude in impact 
parameter space) computed on the basis of DGLAP. We 
obtain an excellent fit of the form a(x + x~y , where x + 
and x~ are the momentum fractions of the "first" ladder 
partons on respectively projectile and target side (which 
initiate the parton evolutions). The parameters a and j3 
depend on the cms energy y/s of the hadron-hadron col- 
lision. To mimic the reduction of the increase of the ex- 
pressions a(x + x~ Y with energy, we simply replace them 

by 

a{x+f +ep (x-f +ET , (I) 

where the values of the positive numbers £p/x will in- 
crease with the nuclear mass number and log s. 

The inelastic rescatterings (ladder splittings, looking 
from inside to outside) amount to providing several lad- 
ders close to the projectile (or target) side, which are 
close to each other in space. They cannot be considered 
as independent color fields (strings), we should rather 
think of a common color field built from several partons 
ladders. We treat this object via an enhancement of rem- 
nant excitations. In fact, the picture described so far is 
not yet complete, since we just considered two interact- 
ing partons, one from the projectile and one from the 
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target. Also the remnants themselves contribute to par- 
ticle production, but mainly in the fragmentation region. 
For more details see Q. 



B. Factorization and Multiple Scattering 

An inclusive cross section is one of the simplest quanti- 
ties to characterize particle production. As discussed ear- 
lier, inclusive cross section are particularly simple, quan- 
tum interference helps to provide simple formulas referred 
to a "factorization". If we want to study high multiplicity 
events, we have to go beyond the inclusive treatment. 

To formulate a consistent multiple scattering theory is 
difficult. A possible solution is Gribov's Pomeron calcu- 
lus, which can be adapted to our language by identifying 
Pomeron and parton ladder. Multiple scattering means 
that one has contributions with several parton ladders 
in parallel. This formulation is equivalent to using the 
eikonal formula to obtain the total cross section from the 
knowledge of the inclusive one. 

We indicated several years ago inconsistencies in this 
approach, proposing an "energy conserving multiple scat- 
tering treatment" The main idea is simple: in case 
of multiple scattering, when it comes to calculating par- 
tial cross sections for double, triple ... scattering, one has 
to explicitly care about the fact that the total energy has 
to be shared among the individual elementary interac- 
tions. In other words, the partons ladders which happen 
to be parallel to each other share the collision energy, 
see fig. [3] A consistent quantum mechanical formulation 
of these ideas requires not only the consideration of the 
usual (open) parton ladders, discussed so far, but also of 




Figure 3: (Color online) Multiple scattering with energy shar- 
ing. 



closed ladders, representing elastic scattering. These are 
the same closed ladders which we introduced earlier in 
connection with elastic rescatterings. The closed ladders 
do not contribute to particle production, but they are 
crucial since they affect substantially the calculations of 
partial cross sections. Actually, the closed ladders sim- 
ply lead to large numbers of interfering contributions for 
the same final state, all of which have to be summed up 
to obtain the corresponding partial cross sections. It is 
a unique feature of our approach to consider explicitly 
energy-momentum sharing at this level (the "E" in the 
name EPOS). For more details see (TT| . 

III. HYDRODYNAMIC EVOLUTION 

A. Parton-Ladders, Flux- Tubes, 
Energy-Momentum Tensor 

In case of high multiplicity pp scattering, we apply ex- 
actly the same procedure as we did for AuAu collisions at 
RHIC, as explained in detail in Q, and shortly reviewed 
in the following. We will identify parton ladders with 
elementary flux tubes, the latter ones treated as classi- 
cal strings. We use the simplest possible string: a two- 
dimensional surfaces X(a,/3) in 3+1 dimensional space- 
time, with piecewise constant initial conditions, referred 
to as kinky strings. In fig. 0Ja), we sketch the space 
components of this object: the string in Et space is a 
mainly longitudinal object (here parallel to the z-axis) 
but due to the kinks (associated to transversely moving 
gluons) there are string pieces moving transversely (in in- 
direction in the picture). But despite these kinks, most 
of the string carries only little transverse momentum! 

In case of elementary reactions like electron-positron 
annihilation or proton proton scattering (at moderately 
relativistic energies), hadron production is realized via 
string breaking, such that string fragments are identified 
with hadrons. When it comes to heavy ion collisions or 
very high energy proton-proton scattering, the procedure 
has to be modified, since the density of strings will be so 
high that they cannot possibly decay independently. For 
technical reasons, we split each string into a sequence of 
string segments, at a given proper-time To, correspond- 
ing to widths 5a and 6(3 in the string parameter space 
(see fig. BJb)). One distinguishes between string seg- 
ments in dense areas (more than some critical density po 
of segments per unit volume), from those in low density 
areas. The high density areas are referred to as core, 
the low density areas as corona [2l|. String segments 
with large transverse momentum (close to a kink) are 
excluded from the core. Based on the four-momenta of 
infinitesimal string segments, 

{ dX(a,l3) s , dX(a,p) sa \ 
Sp= \^^ Sa + ^a— S P}> (2) 
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(a) 



(b) 




Figure 4: (Color online) (a) Flux tube with transverse kink 
in 1R 3 space. The kink leads to transversely moving string 
regions (transverse arrow), (b) String segment at given proper 
time. 



with g being a Gaussian smoothing kernel, one com- 
putes the energy-momentum tensor and conserved cur- 
rents. The corresponding energy density e(ro,x) and the 
flow velocity v(to,x) serve as initial conditions for the 
subsequent hydrodynamic evolutions. 

In fig. [5l we show as an example the energy density 
at ro = 0.6fm/c for a high multiplicity pp collision at 
900 GeV, where high multiplicity here refers to a plateau 



height dn/drj of 12.9, which is more than 3 times the av- 
erage. We see a maximum energy density of about 50 
GeV/fm 3 , which indeed correspond to the energy densi- 
ties observed in central gold-gold collisions at 200 GeV. 
Even more, comparing with the spiky single event results 
for gold-gold in [8j , our pp distribution correspond to one 
(of many) spikes in gold-gold at 200 GeV, which means 
a hydrodynamic treatment for pp is as good (or bad) as 
for gold-gold at 200GeV. 



B. Collective expansion 

Having fixed the initial conditions, matter evolves ac- 
cording to the equations of ideal hydrodynamics, namely 
the local energy-momentum conservation 

9 M T^ = 0, T>™ = (e+p) u^u v - p g^ , (3) 

and the conservation of net charges, 

dN£ = 0, JV£ = n fc ii", , (4) 

with k = B, S, Q, where B, S, and Q refer to respec- 
tively baryon number, strangeness, and electric charge, 
and with u being the four- velocity of the local rest frame. 
Solving the equations, as discussed in the appendix of 
[H, provides the evolution of the space-time dependence 
of the macroscopic quantities energy density s(x), collec- 
tive flow velocity v(x), and the net flavor densities n,k{x). 
Here, the crucial ingredient is the equation of state, which 
closes the set of equations by providing the e-dependence 
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Figure 5: (Color online) Initial energy density in a high multi- 
plicity pp collision (dn/dr] = 12.9) at 900 GeV, at a space-time 
rapidity r] s = 0. 
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Figure 6: (Color online) Energy density versus temperature, 
for our equation-of-state X3F (full line), compared to lattice 
data [22l | (points), and some other EoS choices, see @|. The 
thin vertical line indicates the "hadronization temperature" 
Th, i.e. end of the thermal phase, when "matter" is trans- 
formed into hadrons. 
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Figure 7: (Color online) Energy density (upper panel) and radial flow velocity (lower panel) for a high multiplicity pp collision 
(dn/dr] = 12.9) at 900 GeV, at proper times r = 1.3 fm/c (left) and r = 1.9 fm/c (right), at a space-time rapidity r) a = 0. 



of the pressure p. As discussed in we use an equa- 
tion of state compatible with lattice gauge simulations, 
see fig. |U 

Starting from the flux-tube initial condition, the sys- 
tem expands very rapidly. It hadronizes in the cross-over 
region, where here "hadronization" is meant to be the end 
of the completely thermal phase: matter is transformed 
into hadrons. We stop the hydrodynamical evolution at 
this point, but particles are not yet free. Our favorite 
hadronization temperature is 166 MeV, shown as the thin 
vertical line in fig. [51 which is indeed right in the tran- 
sition region, where the energy density varies strongly 
with temperature. At this point we employ statistical 
hadronization, which should be understood as hadroniza- 
tion of the quark-gluon plasma state into a hadronic sys- 



tem, at an early stage, not the decay of a resonance gas 
in equilibrium. 

After this hadronization -although no longer thermal- 
the system still interacts via hadronic scatterings. The 
particles at their hadronization positions (on the corre- 
sponding hypersurface) are fed into the hadronic cascade 
model UrQMD [23|, [24j , performing hadronic interactions 
until the system is so dilute that no interactions occur 
any more. The "final" freeze out position of the particles 
is the last interaction point of the cascade process, or the 
hydro hadronization position, if no hadronic interactions 
occurs. 

In fig. [JJ we show the hydrodynamic evolution of the 
event corresponding to the initial energy density of fig. 
which can be considered as a typical example, with simi- 
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lar observations being true for randomly chosen events of 
this multiplicity (dn/dr] — 12.9). We see that the system 
evolves immediately also transversely, the energy density 
drops very quickly. A very large transverse flow develops, 
typically around 70 % of the velocity of light. This will 
have measurable consequences. 



IV. ELEMENTARY DISTRIBUTIONS 

We first check some elementary distributions. We use 
the EPOS 2.05 version, which has been optimized for 
heavy ion scattering at RHIC, the same one as used in 
Q . We could certainly improve the results by doing some 
"tuning" taking into account the new LHC results, but the 
purpose of this paper is more to show what we get from 
a straight application of the "heavy ion model", here ap- 
plied to pp at LHC. We only consider 900 GeV, for higher 
energies some reconsideration of our screening procedures 
will be necessary (work in progress). As usual we work 
with the event-by-event mode, and hydrodynamics is only 
employed for high density areas (core-corona separation). 

In the following we will compare three different scenar- 
ios: 

full : the full calculations, including hydro evolution and 
hadronic cascade; 

no case : calculation without hadronic cascade; 

base ; calculation without hydro and without cascade. 

We will compare the corresponding calculations with ex- 
perimental data, for pp scattering at 900 GeV. 
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Figure 8: (Color online) Pseudorapidity distributions in pp 
scattering at 900 GeV, compared to data (points). We show 
the full calculation (solid line), a calculation without hadronic 
cascade (dashed), and a calculation without hydro and with- 
out cascade (dotted). 
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Figure 9: (Color online) Transverse momentum distributions 
in pp scattering at 900 GeV, for minimum bias events (up- 
per panel) and high muliplicity events (n = 22, lower panel), 
compared to data (points). We show the full calculations (solid 
lines), a calculation without hadronic cascade (dashed), and 
a calculation without hydro and without cascade (dotted). 



In fig. [5J we show pseudorapidity distributions of 
charged particles, compared to data from CMS [25[ and 
ALICE [26|, [2tJ ■ The three scenarios do not differ very 
much, and agree roughly with the data. 

We then investigate transverse momentum distribu- 
tions. For minimum bias events, there is again little dif- 
ference for the three scenarios (all of them reproduce the 
data within 20%), as seen in the upper panel of fig. [5] 
The situation changes drastically, when we consider high 
multiplicity events, see the lower panel of fig. [9j Here 
the base calculation (without hydro) underestimates the 
data by a factor of three, whereas the full calculation gets 
close to the data. This is a very typical behavior of col- 
lective flow: the distributions get harder at intermediate 
values of p t (around 1-4 GeV/c). 

In fig. [TU1 we plot the mean transverse momentum as 
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Figure 10: (Color online) Mean transverse momentum as a 
function of the charged multiplicity in pp scattering at 900 
GeV, compared to data (points). We show the full calcula- 
tion (solid line), and a calculation without hydro and without 
cascade (dotted). 



a function of the charged multiplicity, compared to data 
from ALICE [13] • The increase of the mean p t with mul- 
tiplicity is in our approach related to collective flow: with 
increasing multiplicity one gets higher initial energy den- 
sities, and more collective flow can develop. The data are 
therefore compatible with our flow picture, but for a real 
proof one needs at least in addition the mean p t behavior 
of heavier particles (protons, lambdas, or even heavier), 
since the effect gets bigger with increasing mass. 



as expected a steeply falling distribution as a function of 
r. In the two cases involving a hydrodynamical evolution, 
particle production is significantly delayed, even more in 
the case of the full calculation, with hadronic cascade. 
The bump in the two latter scenarios is due to particles 
being produced from the fluid, the small pt contribution 
is due to corona particles. 

This particular space-time behavior of the hydrody- 
namic expansions should clearly affect Bose-Einstein cor- 
relations - what we are going to investigate in the fol- 
lowing. There is a long history of so-called femtoscopic 
methods [28 32], where the study of two-particle corre- 
lations provides information about the source function 
S(P, r'), being the probability of emitting a pair with to- 
tal momentum P and relative distance r'. Under certain 
assumptions, the source function is related to the mea- 
surable two-particle correlation function CF(P,q) as 



CF(P,q) = J dVS(P,r')|*(q',r') 



(6) 



with q being the relative momentum, and where ^ is 
the outgoing two-particle wave function, with q' and r' 
being relative momentum and distance in the pair center- 
of-mass system. The source function S can be obtained 
from our simulations, concerning the pair wave function, 
we follow [HI, some details are given in Q. 

Here, we investigate tt + - tt + correlations. We evaluate 
eq. (J5]), with Bose-Einstein (BE) quantum statistics in- 
cluded, but no Coulomb corrections. Weak decays are not 
carried out. In figs. [T2lll3ll^l we show the results for dif- 
ferent k T intervals defined as (in MeV): KT1= [100, 250], 



V. BOSE-EINSTEIN CORRELATIONS 

The space-time evolution of the "full" hydrodynamic 
approach will be completely different compared to the 
"base" approach, where particles are directly produced 
from breaking strings, as can be seen from fig. 1111 where 
we plot the distribution of formation points of ir + as a 
function of the radial distance 

r = \J x 2 + y 2 (5) 

(in the pp center of mass system (cms)). Only parti- 
cles with space-time rapidities around zero are consid- 
ered. We compare aain the three scenarios "full" (full 
calculation - flux-tube initial conditions, hydro, hadronic 
cascade), "no case" (without hadronic cascade, only flux- 
tube initial conditions and hydro, hadronization as usual 
at 166 MeV), and "base" (without hydro and without cas- 
cade, just flux-tube approach with string decay). 

All calculations in this section refer to high multiplicity 
events in pp scattering at 900 GeV, with a mean dn/dr](0) 
equal to 12.9. The "base calculation" (dotted line) gives 




0.5 1 1.5 2 2.5 3 
radial distance r (fm) 



Figure 11: (Color online) The distribution of formation points 
of 7r + as a function of the radial distance in a high multiplicity 
event from pp scattering at 900 GeV, for the following sce- 
narios: the full calculation (solid line), a calculation without 
hadronic cascade (dashed), and a calculation without hydro 
and without cascade (dotted). 
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KT3= [400,550], KT5= [700, 1000], where k T of the pair 
is defined as 



k T = g (lft(P ionl ) 



p t (pion2)|) . 



(7) 



We compare the three different scenarios: "full calcula- 
tion" (solid line), "calculation without hadronic cascade" 
(dashed), and "calculation without hydro and without 
cascade" (dotted), and data from ALICE [34| . The data 
are actually not Coulomb corrected, because the effect 
is estimated to be small compared to the statistical er- 
rors. We consider here the high multiplicity class, with 
dn/d-q(0) = 11.2, close to the value of 12.9 from our sim- 
ulated high multiplicity events. We compare with the 
real data (not polluted with simulations) , normalized via 
mixed events, and we do the same with our simulations. 
Despite the limited statistics, in particular at large kx, 
we see very clearly that the "full" scenario, including hy- 
dro evolution and hadronic cascade, seems to fit the data 
much better then the two other ones. Usually people 
like to extract radii from these distributions, so when we 
make a fit of the form 



CF — 1 = A cxp(-i?|q|), 



(8) 



in the |q| range from 0.05 to 0.70. We obtain the radii 
given in the figure. So the radii are very different, vary- 
ing from 0.69 fm (base approach) to 1.80 fm (full model), 
which is understandable from fig. [TTJ We prefer an ex- 
ponential fit rather than a Gaussian, simply because the 
former one works, the latter one does not. We do not 
want to give a precise meaning to R, it simply character- 
izes the distribution. 

Normalizing by mixed events is something one can eas- 
ily do experimentally (this is why we compare with these 
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Figure 12: (Color online) The correlation functions CF for 
7r + - 7r + pairs as obtained from our simulations, for the three 
different scenarios, for kr bin KT1, compared to data (points). 
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Figure 13: (Color online) Same as fig. 1121 but kr range KT3. 



data), but it is clear that one has still unwanted cor- 
relations, like those due to energy-momentum conserva- 
tion, which is not an issue in mixed events. Doing sim- 
ulations, life is easier. We can take simulations without 
Bose-Einstein correlations as base line, rather than mixed 
events. This is referred to as "real / bare" normalization 
(to be distinguished from the "real / mixed" case dis- 
cussed earlier). The corresponding results are show in 
fig. HSl the solid line (full calculation) is now completely 
horizontal away from the peak region, the radius from 
the exponential fit is 2.10 fm instead of 1.80 fm for the 
"mixed" normalization. For the other kx regions, the sit- 
uation is similar, the final results for all three kx regions 
for the full calculation is shown in fig. [151 together with 
the the radii from the exponential fit: they are almost 
identical, around 2 fm. 

We get to the same conclusion as outlined in [34j : the 
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Figure 14: (Color online) Same as fig. 1121 but kr range KT5. 
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Figure 15: (Color online) Same as fig. 1121 but normalization 
via a simulation w/o BE correlation ("bare"). 



Figure 17: (Color online) As fig. 1161 but calculation without 
hydro and without cascade. 



radii are kx independent, contrary to what has been ob- 
served in AuAu scattering. 

How can it be that our hydrodynamic scenario gives 
a strong kx dependence in AuAu, but not in pp? To 
answer this question, we compute the "true" correlation 
function (real / bare normalization) for the calculation 
without hydro and without cascade (just string decay). 
The results are shown in fig. 1171 Surprisingly, here we get 
a strong kx dependence of the radii, but the "wrong" way: 
we have 0.64 for KT1 and 1.63 fm for KT5 ! Actually 
such a behavior is quite normal, as seen from fig. [TBI 
the distribution is broader for high p t particles, because 
high p t resonances live longer and can move further out 
before decaying. This effect is in principle also present 
in AuAu scattering, but it is much more visible for the 



small pp system. So in pp we have two competing effects: 

• radii increase with kx , due to the bigger size of the 
source of the high p t particles compared to the low 
Pt ones, 

• radii decrease with kx , as in AuAu (see Q ) , in case 
collective flow, due to the p-x correlation. 

As seen in fig. [lj|l this p — x correlation exists indeed for 
the case of hydrodynamic evolutions, and is much smaller 
in the basic scenario. So in the hydro scenarios, the two 
competing effects roughly cancel, the radii are kx inde- 
pendent. To really see the x — p correlation, one need 
to "divide out" the trivial kx dependence due to the pt 




q 

Figure 18: (Color online) The distribution of formation points 
Figure 16: (Color online) Correlation function, normalized by of particles as a function of the radial distance, for the scenario 
using a simulation w/o BE correlation ("bare"), for three k T "without hydro and without cascade", for high p t and low p t 
intervals. particles. 
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Figure 19: (Color online) Momentum - space correlation, for 
the different scenarios. 



dependence of the single particle source sizes, which we 
do by considering the kx dependence of R/R\, as , with 
the reference radius i?bas referring to the base scenario 
(without hydro, without cascade), see fig. [20] The ra- 
tio R/Rhas decreases with kx as a manifestation of the 
x — p correlation, as a consequence of the hydrodynamic 
expansion. 

An alternative way of getting out unwanted correla- 
tions would be the consideration of double ratios like 

CF(full scenario with BE) /CF(full scenario w/o BE) 
CF(base scenario with BE) /CF(base scenario w/o BE) ' 

(9) 

where basic scenario refers to the calculation without hy- 
dro and without hadronic cascade. 
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k T (GeV/c) 

Figure 20: (Color online) The kx dependence of R/Rb aB . The 
ratio decreases significantly with kr, a clear "flow signal". 



VI. SUMMARY 



After having introduced recently a sophisticated ap- 
proach of hydrodynamic expansion based on flux-tube 
initial conditions for AuAu collisions at RHIC, we now 
employ exactly the same picture to pp scattering at 900 
GeV, which is in particular justified for high multiplicity 
events. A very interesting application are Bose-Einstein 
correlations. We have shown that as in heavy ion scat- 
tering the hydrodynamic expansion leads to momentum 
- space correlations, which clearly affect the correlation 
functions. To see the signal is non-trivial due to the fact 
that in addition to the x — p correlations (which leads to 
decreasing radii with kx, there is a second effect which 
works the other way round: the single particle source size 
is p t dependent, which is an important effect in pp, not 
so in heavy ion scattering. In this sense we can interpret 
the kx independence of the radii as a real flow effect. 
Our simulation does not only reproduce the kx indepen- 
dence, but also the whole correlation functions, which is 
not at all reproduced from the "base scenario" without 
hydro and without cascade. So the correlation data pro- 
vide a very strong evidence for a collective hydrodynamic 
expansion in pp scattering at the LHC. 
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